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ABSTRACT 

When the ion mean free path much exceeds the Larmor radius in a plasma, 
the viscous stress tensor is altered dramatically, and depends only upon quanti- 
ties measured along the field lines. This regime corresponds to typical interstellar 
medium conditions in galaxies and protogalaxies, even if the magnetic field is ex- 
tremely weak, with a negligible Lorentz force on all scales of interest. In this 
work, the only role of the magnetic field is to channel angular momentum trans- 
port along its lines of force. We show that differential rotation in such a gas is 
highly unstable, with a maximum growth rate exceeding that of the magnetoro- 
tational instability. The regime of interest has been treated previously by plasma 
kinetic methods. Where there is overlap, our work appears to be in agreement 
with the kinetic results. The nonlinear outcome of this instability is likely to be 
a turbulent process, significantly augmenting the magnetorotational instability, 
and important to the initial phases of the amplification of small galactic magnetic 
fields. 
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1. Introduction 

The magnetorotational instability, or MRI, has become central to our understanding 
of turbulent angular momentum transport in accretion disks (e.g. Balbus 2003). The in- 
stability is important even when (indeed, especially when) the magnetic energy density is 
small compared with the thermal energy density. It is this feature, the fact that the stability 
of the gas is hypersensitive to the presence of subthermal magnetic fields of any geometry, 
that endows the MRI with its special significance. In fact, the MRI is merely one mani- 
festation of much more general behavior. By imparting new degrees of freedom to a fluid, 
magnetic fields allow free energy gradients to become sources of instability, with important 
consequences for a variety of astrophysical systems: accretion disks become turbulent when 
the angular velocity (not angular momentum) decreases outward, and stratified dilute gases 
are destabilized when the temperature (not entropy) decreases upwards (Balbus 2001). 

The thermal destabilization caused by magnetic field is especially noteworthy because 
it is engendered by what seems at first a purely dissipative process: thermal diffusivity. 
It is the extreme anisotropy of the conductivity tensor parallel and perpendicular to the 
field lines that lies at the heart of the instability. By channeling heat only along lines of 
magnetic force, small-amplitude ripples along initially isothermal field lines grow into large 
fluid displacements parallel to the temperature gradient. This is not, as is sometimes thought, 
analogous to classical double-diffusive instabilities, such as ocean layer "salt-fingering" 1 . 
Rather, it depends wholly upon the properties of anisotropic conductivity. 

In this paper we show that the anisotropy introduced into the viscous stress tensor by 
a weak magnetic field sharply destabilizes dilute astrophysical disks, even without Lorentz 
forces appearing in the fluid equations. By "dilute," we mean the limit in which the ion 
Larmor radius is small compared to a mean free path, which in turn is small compared with 
the characteristic macroscopic length scales of the disk. We shall refer to the resulting viscous 
stress tensor in this limit as Braginskii viscosity (Braginskii 1965). The Braginskii limit is 
appropriate for interstellar and galactic disks (especially protogalactic disks, cf. Malyshkin 
& Kulsrud 2002). A remarkable property of the instability is that in the absence of Lorentz 
forces, when the viscous diffusivity much exceeds the resistive diffusivity, rapid growth times 
are associated with arbitrarily high wave numbers. Except for isolated field geometries (e.g. 
precisely azimuthal), there is no formal high wavenumber dissipation of the linear magneto- 
viscous instability. By way of contrast, for the ultra-weak magnetic fields considered here, 
the wavelength of maximum growth in the classical MRI would be strongly damped by an 
ordinary isotropic viscosity. (Note that in real systems, at sufficiently high wave numbers, 



1 See Mcnou, Balbus, & Spruit (2004) for a true salt-fingering analogy involving the MRI. 
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Lorentz forces will ultimately stabilize.) 

The magnetic stability of a dilute plasma may also be studied using a plasma kinetic 
approach. This offers rigor in a problem where it is clearly of some benefit, but at a price 
of greater mathematical complexity. Quataert, Dorland, & Hammett (2002) analyzed the 
MRI in the collisionless regime, and noted that the character of the instability changes when 
the pressure tensor becomes significantly anisotropic, with growth rates in excess of the 
classical MRI maximum. In the current work, we have in essence abstracted the anisotropic 
components of the pressure, and labeled them as a viscosity tensor. More recently, Sharma, 
Hammett & Quataert (2003) reanalyzed the kinetic problem using a Krook collision operator, 
and showed that the collisional and collisionless behavior are on the same branch of the 
dispersion relation. In the short mean free path limit, the pressure tensor associated with 
the equations used by Sharma et al. 2003 reduces to a scalar pressure plus the dominant 
parallel components of the Braginskii viscosity (Snyder, Hammett, & Dorland 1997). The 
current work, a purely gasdynamical treatment of the problem, affords relative mathematical 
simplicity and a physically transparent interpretation. 

In either its kinetic or gasdynamical guise, this vigorous instability may be important 
in the early stages of magnetic field amplification in disk galaxies, when densities are low, 
temperatures are relatively high, and the field is likely to be very weak. Quataert et al. 
(2002) and Sharma et al. (2003) discuss applications to black hole accretion flows of low 
radiative efficiency 

An outline of the paper is as follows. Section 2 is a discussion of the regime of the 
applicability of this work and a presentation of the magnetic viscosity formalism. Section 3 
is the heart of the paper, formulating and solving the problem, and checking the validity of 
the results. Section 4 is a discussion of the applicability of our results to galactic magnetism, 
and section 5 summarizes our findings. 

2. Preliminaries 

2.1. Gasdynamical Description of the Instability 

We begin with a physical description of the instability. The process is very simple. In 
the leftmost diagram of figure 1 (a), we show three representative azimuthal field lines, with 
increasing azimuth toward the top. The angular velocity gradient Vf2 is radial, pointing 
from right to left, so that smaller angular velocities lie to the right. The angular momentum 
gradient, by contrast, runs from left to right. Initially, there is no viscous transport of 
angular momentum, because the angular velocity gradient and field lines are orthogonal. 
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The middle portion of figure 1 (b) shows the field lines and two fiducial fluid elements 
(dark circular disks) after a perturbation has created a small oscillatory radial component. 
Now there is a component of Vf2 along the field lines, and viscous angular momentum 
transport, represented by the thin arrows, is able to commence. The rightward fluid element 
acquires angular momentum, while the leftward element loses angular momentum. 

In the rightmost portion of the figure (c), the loss of angular momentum causes the 
element on the left to fall further leftward, towards orbits of smaller angular momentum. 
The element on the right, on the other hand, acquires angular momentum, and moves toward 
the higher angular momentum orbits that lie further to the right. The critical point is that 
these dynamical motions stretch the field lines yet more radially, allowing yet more angular 
momentum to flow between the elements, and the process runs away. 

The instability blends different aspects of previously known instabilities. It works on 
the same dynamical principles as the MRI, tethering fluid elements and transporting angular 
momentum between them, but there is no magnetic tension here. Similarly, it works using 
the same gradient transport process seen in the magneto-thermal instability, channeling 
diffusive transport along magnetic field lines while simultaneously realigning them, but with 
angular momentum rather than heat. The result of all this, we shall see, is an extremely 
powerful shear instability. 



2.2. Dilute Plasma Limit 

The presence of a magnetic field alters the form of the viscosity when 

e = {u ci r ci )~ l < 1, (1) 

where u) ci is the ion cyclotron frequency and r ci is the ion-ion mean collision time. (This is 
tantamount to having the mean free path much exceed the gyroradius.) Under laboratory 
conditions, this is normally considered to be the large field strength limit, but it is very 
much the norm for interstellar plasmas, even when the field is very weak. Indeed, reference 
to Spitzer (1962) gives for a hydrogenic plasma 

1.09 x 10 5 \ Tl /2 B^ G 



( 2) 

where n is the proton density in cm~ 3 , T4 the kinetic temperature in units of 10 4 K, B^q is 
the magnetic field in microgauss, and In A is the Coulomb logarithm, ~ 20. With n < 1 and 
T 4 > 1, very weak fields are clearly accommodated by the asymptotic regime e < 1. The 
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ion mean free path X m f P and gyroradius r g are respectively 

\ mfp = 5.2 x 10 11 (J^ (3) cm, r g = 1.3 x 10 8 ] cm. (3) 



2.3. Ideal MHD Limit 

Throughout this work, we ignore the effects of finite resistivity. The ratio of the viscosity 
to resistivity is known as the magnetic Prandtl number V, and may be adapted from Balbus 
& Hawley (1998): 

T V /6.5 x 10 10 \ / 20 



r ={w) I^^JUaJ • < 4 > 

We may easily restrict the calculations to T > 10 4 K and n <C 10 10 cm -3 , so that V >> 1, 
and resistivity will be ignored. Since these questions can be subtle however, we return to 
this point ex post facto in §3.5. 



2.4. Magnetic Viscous Stress Tensor 

The theory of viscous transport in magnetized plasmas is presented by Braginskii (1965). 
The usual isotropic collisional viscous stress tensor can be written 

o-ij = -rjWij, (5) 

where 77 is the dynamical viscosity coefficient, and 

dvj dvj 2 . , 

This form applies to a set of Cartesian axes, k) being an even permutation of (X, Y, Z). 
As usual, Sij denotes the Kronecker delta function. We note that the stress is traceless. In 
the paper, we work exclusively in the Boussinesq limit, and shall set V-v = in the above. 
In §3.4, it is shown that the Boussinesq limit is justified when the Reynolds number is large. 

In the presence of a restricting magnetic field, the only component of that remains 
unaffected is the momentum flux along the magnetic line of force due to the gradient along 
the field line. Define a local Cartesian coordinate system (the "field frame") (X b ,Y b , Z b ), 
chosen with the magnetic field lying along the Z b axis. Then 

<?z b z b = ^bibjo-ij, (7) 
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where the h are components of the unit magnetic field vector in an arbitrary Cartesian 
frame. Braginskii (1965) shows that all other components of the magnetized viscous stress 
tensor are smaller than uz b z b by a factor of order e or e 2 , and are therefore ignored in this 
calculation. The important exception are the two other diagonal stress components, which 
to leading order in e are identical with one another. Since the stress must always be traceless, 
we have 

°X b X b = °~Y b Y b = --CZ b Z b - (8) 

In vector notation, the Z^Zf, component of the stress is then given by 



<Jz b z b = -277 l(b-V)v] -b. (9) 

To find the components of the magnetized viscous stress tensor in any other locally Cartesian 
frame, the transformation law may be written 

a H = J^ii-ib) tf-jb) °i b j b , (10) 

ib,jb 

where once again the b subscript denotes the magnetic field frame and bold face quantities are 
unit vectors of the indicated component. Using equation (8) for the nonvanishing diagonal 
stress tensor components in the field frame, we find 



- <5 Zh Z h 



[i-Z b )(j-Z b ) - ±(i.Y b )(j.Y b ) - l -{i-X b ){j-X b ) 



;n) 



Equations (9) and (11) allow one to determine the magnetic viscous stress in a frame suitable 
for working with fluid variables, given the local field geometry. 



3. Formulation of the Problem 

We consider the stability of a disk in the presence of a magnetic field, but with a field 
so weak that all dynamical magnetic forces are negligibly small. In contrast to the MRI, we 
assume this "magneto-anemic" condition is true not only for the equilibrium disk, but even 
for small wavelength, WKB perturbations. The sole effect of the magnetic field is to restrict 
viscous transport in accordance with prescription of §2. 

The fundamental fluid equations used here are mass conservation 

^ + V-M = 0, (12) 

the equation of motion, 

( ,(| + „. V )«=-VP-pV»-g, (13) 
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and the induction equation of ideal MHD, 



dB 



Vx(t)XB). 



(14) 



dt 



For reasons that will become clear, an internal energy equation is not needed at this stage. 

The equilibrium state is a differentially rotating disk, and we work in the usual cylindri- 
cal coordinate system, R, <p, Z. The angular velocity is Q(R), and we shall restrict ourselves 
to a local analysis at the midplane. Thus, we may ignore buoyant forces, which depend upon 
gradients in pressure and entropy. 

In the equilibrium state, it is assumed that cr^- = 0. Indeed, the point of this calculation 
is to show that any development of viscous transport along a field line is highly unstable. 
The initial magnetic field lines are wrapped around cylinders, and are unaffected by the 
shear. 

We consider next small departures from the equilibrium flow. Linearly perturbed quan- 
tities are denoted by dv,5<Jij, etc. We work in the local WKB limit, with the space-time 
dependence of all perturbed quantities given by exp^t + ik-r). Thus, 7 is a growth or decay 
rate if it is real, and an angular frequency if it is imaginary. 

The wavenumber k as well as the assumed constancy of 7 require some further ex- 
planation. The wavenumber has radial, azimuthal and vertical components kn,m/ R, k z 
respectively. Since we are working in a local shearing system, the radial wavenumber Ur will 
formally depend on time (Goldreich & Lynden-Bell 1965; Balbus & Hawley 1992): 



where &r(0) is the initial value of k^. For our present purposes however, the time dependence 
of kn will prove irrelevant, as the radial wavenumber disappears early in the analysis. It is 
for this reason that we may also assume a simple exponential time dependence; in general 
the problem is more complex. 

To evaluate 5<Jij, it follows from equation (11) that 



since cr ZbZb vanishes in the unperturbed state. The terms in square brackets may be evaluated 
for the equilibrium field geometry, which we take to be 



k R (t) = fcfl(0) - mt 



(15) 



d\nR 



6** = Sa ZbZb (i-Z b )(j.Z b ) - l -(i.Y b ){j-Y b ) - l -{i-X b ){j-X b ) 



(16) 



b = cos 9 4> + sin 9 Z. 



(17) 
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<f>, Z, and R are unit vectors in the indicated cylindrical directions, and 9 is the angle 
between the magnetic field and the axis. The local magnetic field axes are then 



Z b = b = cos 9 4> + sin 9 Z, 
X b = Rx Z b = - sin 9 <fi + cos 9 Z, 



(18) 
(19) 



and 



Y b = R. (20) 

See fig. 2. 

The nonvanishing elements of Sa^ may now be calculated. The diagonal elements are 



5&rr = --5&z b z b , 5a 4 



cos 2 9 - 



sin 2 9 



$°~z b z b , Sa Z z = sin 2 9 - 



cos 2 9 



5a 



and the off-diagonal elements are 

3 

5cr<pz = Suzcp = - cos 9 sin 9 5a ZbZb . 
We may evaluate 5az b z b from its vector-invariant form (9), 

Sa ZbZb = — 277 ([(Sb-V)v] -b + [(b-V)Sv] -b + [(b-V)v] -6b) 
Using v = RVLcf) and equation (18), we find 

dVt 



Sa Zb z b = -2r] 



5b 



R 



d\nR 



cos9 + i(k-b)(b-Sv) 



Z b Z b j 

(21) 
(22) 

(23) 
(24) 



3.1. Linearized Equations 

The linearized equations are written in a coordinate system that is shearing with the 
unperturbed flow. The only effects this has on the equations of motion are that the time 
derivative must be Lagrangian, 

* s a? + V (25 > 

and the radial partial derivative is replaced by the kn(t) wavenumber in equation (15) (Balbus 
& Hawley 1992). The dynamical equations are 



k-dv = 0, 



(26) 
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Fig. 2. — Relative orientation of cylindrical R, 0, Z and magnetic axes. Magnetic 

field lies along and R and Y& axes coincide. 
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^-msv^ + ^(sp- 5 -^ 

dt P V 2 



(27) 



d 8vs k 2 . m 



SP ( cos 2 0-^1 <)a 



ZbZi, 



.3k 



+ i—^sm9 cos 9 5a ZbZb = 0, (28) 



2p 



dbv z ikz_ 
dt p 



S P 4- |sm 2 - — ) 



. 3m 



sin# cos#<5a ZbZb = 0. (29) 



In equation (28), we have introduced the standard notation k for the epicyclic frequency, 
which characterizes the response of fluid element displacements in the R<p plane in a non- 
magnetized disk. In terms of the angular velocity fl(R), 



k = 4tt 2 + 



dtt 2 
dlnR' 



(30) 



Notice that only the perturbed radial component of the magnetic field unit vector Sbn 
is required in 5az b z b - This is easily obtained from the radial induction equation: 



dSh 
~dt 



i(k-b) Svr. 



(31) 



When the time dependence of the perturbations is a simple exponential, this immediately 
gives 

i(h..h\ 

(32) 



oo R = 5v R , 



and 



8&z b z b = -2r)i(k-b) cos6> 



7 

5vr dfl 
7 dlnR 



+ Sv d 



— 2r/i(k-b) s\n98vz- 



(33) 



3.2. Axisymmetric Disturbances 

If there are both vertical and azimuthal field components present, the problem admits 
a relatively simple axisymmetric solution with k having only a vertical component. This is 
also likely to be the most unstable mode, and therefore of greatest astrophysical interest. 



With k = kzz, mass conservation immediately gives 8vz = 0. The remaining dynamical 
equations are 

d 8vr 

2JWu^ = 0, (34) 



dt 
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SP (sin 2 ^1 :(). (36) 



Equation (36) is one of vertical hydrostatic equilibrium. In an adiabatic gas, this constraint 
must be consistent with the Boussinesq approximation, a point we revisit in §3.4. 

The problem decouples, and we are simply left with the two equations (34) and (35) 
for 5vr and Sv^. We may seek solutions with an exponential time dependence of the form 
exp(7t). Then 

^ = -2vik z sin 9 cos ( J^Ltji + Sv ) , (37 ) 
p \ a In R 7 J 

where we have introduced the kinematic viscosity, 

v = ^. (38) 
P 

Substituting this into equations (34) and (35) yields the dispersion relation 

dtt 2 

7 3 + (3uk 2 z sin 2 9 cos 2 #h 2 + k 2 7 + Zvk 2 z sin 2 cos 2 9 — — = 0. (39) 

d In i? 

This may be written 

3^f sin 2 9 cos 2 = - ^It/f? o > (40) 

from which it is clear there is an unstable branch (7 > 0) of the dispersion relation if 
7 + dQ 2 /dR < 0. The maximum growth rate is 

2 dn 2 

lmax d\nR' [ ] 

It is noteworthy that 7 may exceed the Oort A value (l/2)\dQ/dlnR\, a result emphasized 
by Quataert et al. (2002) for a collisionless gas, who also recovered equation (41) in the 
appropriate limit. We note here the interesting fact that the growth rate given in (41) will 
always exceed the Oort A value in any disk that is locally stable by the Rayleigh criterion. 
Equations (34) and (35) may be compared with equations (34) and (35) of Quataert et al. 
(2002). Note in particular the identification of the 4>Z component of the viscous stress tensor 
with the term proportional to Sp\\ — 5p± on the right side of equation (35) of Quataert et 
al. (2002). In both the kinetic and gasdynamical treatment, the azimuthal torque has no 
stabilizing radial counterpart, which is responsible for the the resulting growth rate exceeding 
the classical MRI Oort A value. 

For a typical flat galactic rotation curve, the above gives 7 ma:r = \/2VL. This is an 
enormous rate: linear amplitudes grow a factor of 7.2 x 10 3 in one orbit, 5.2 x 10 7 in two 
orbits. We have of course neglected magnetic tension, which would ultimately be a stabilizing 
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influence, but not before it was a powerful destabilizing influence in the form of the MRI. 
A general treatment including the dynamical effects of the field is deferred to a forthcoming 
paper. But our finding certainly suggests that a kinematic treatment of field amplification 
in galactic disks is at best questionable. For plasma kinetic treatments including Lorentz 
forces, see Quataert et al. (2002) and Sharma et al. (2003). 



3.3. Azimuthal Field Instability 

The above analysis depended upon both axial and azimuthal field components being 
present. It is of interest to isolate the case of an exactly azimuthal field {6 = 0) and examine 
its stability properties. 

We assume that the wave vector is dominated by its axial component; however, in 
contrast to the previous section, we do not set the azimuthal and radial components to zero. 
Rather, we work in the asymptotic limit 

k z ^> k R {0),m/R. (42) 

The dynamical equations of motion in this case are 

^- 2Mtv + ^?(^_^) =0 , (43) 



d5v, 



4> 



. III 



dS »z + *L( ' SP - S -ZMl] = o. (45) 



dt p 

Mass conservation k-dv = implies that 5v z is of order 1/kz relative to the planar velocity 
components, and equation (45) simplifies to 

5P - 5a ZbZb = 0. (46) 

The remaining dynamical equations are 

dSv R 



dt 



2Q5V4, = 0, (47) 



dSvd, .m{3. i 
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Once again, the radial wavenumber Ur has dropped from the problem (as has the dominant 
component kz), and we may look for simple exponential time behavior of the form e 7 *. We 
find 



rn 



5vr dVt 
7 dlnR 



+ Sv d 



(49) 



But this is exactly the problem we solved in the previous section, with the substitution 

kz sin 9 cos 6 <— mj R. (50) 
We obtain exactly the same dispersion relation, 

dQ 2 



Y + 3Yvm 2 /R 2 + k 2 1 + Su(m /R ) J l .. .. 

duiR 

and thus exactly the same maximum growth rate, \dQ 2 /d\nR\. 



0, 



(51) 



3.4. Validity of the Boussinesq Approximation 

The Boussinesq approximation was very important for simplifying the preceding analy- 
ses. We review its validity here. 

The key point is that the perturbed pressure is of order 

5P ~ 5a ZbZb ~ (k-b) V 5v^ (52) 

since 8vr and Sv^ are comparable and the growth rates are of order Q. Assuming that any 
radial entropy gradient is very small, 

bp 5P 5P 



P p pcV 
where c$ is the isothermal sound speed. Hence 

5p u(k-b) 5v<p 



P c s c s 
The Boussinesq approximation requires 



(53) 



(54) 



dhvp ^5p , . 

—1- ~ fi— < k5v, (55) 
dt p 

where k and 5v are respectively characteristic values for the wave number and perturbed 
velocity. With k ~ k-b, this requirement simplifies to 



v 



(4/n) " 1 (56) 

The quantity on the left is 1/Re, the inverse of the Reynolds number. The Boussinesq limit 
is therefore generally appropriate for our problem in the limit of large Reynolds number. 
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3.5. Effect of Resistivity 

A finite resistivity will always be present, directly affecting the dynamics because of its 
role in modifying the constraint of field freezing. The question arises as to whether resistivity 
may generally be ignored at all wavenumbers if the viscous diffusivity v much exceeds the 
resistive diffusivity i] B . (Unlike the transformation of the viscous stress tensor, the parallel 
and transverse resistivities do not differ profoundly [Spitzer 1962], and we shall ignore the 
distinction here.) 

The presence of Ohmic resistance alters equation (32) to 

i(k-b) 



7 + r] B k 2 



5v R , (57) 



where k 2 is the magnitude of the wavenumber. In both of our worked examples, k may be 
taken as kz- 

It is a straightforward exercise to rework the dispersion formula using the above for 6r. 
One obtains for the case of axisymmetric disturbances, 

7 3 + i 2 k 2 z {rj B + 3v sin 2 9 cos 2 9) + 7(/t 2 + 3i] B vk\ sin 2 9 cos 2 G) + 

+k 2 z t] B k 2 + 3v sin 2 9 cos 2 9— — — = 0. (58) 



dlnR 



As an equation for k 2 z this reads 



34(7^ sin 2 9 cos 2 9) + k 2 z ^ B ( 7 2 + k 2 ) + 3v sin 2 9 cos 2 9( 1 2 + dVt 2 jd In R) 

+ 7 ( 7 2 + k 2 ) = 0. (59) 

For 7 > 0, positive solutions for k z exist only if the coefficient of k 2 z is negative, which, in 
the limit v ^> r] B , requires dQ 2 /dR < as before. The quadratic discriminant in the solution 
for k z must also be positive. If we set 

2 dQ 2 dfl 2 . , 

7 + d^R = e ^ (60) 

and assume that e < 1, the condition for the positivity of the discriminant is easily worked 
out. We find 



e2> 64,, 



dlnfl 2 



dlnR 



(61) 



3z/sin 2 (2fl) 

Provided that 9 is not too small, this is consistent with our small e assumption, and shows 
that a small resistivity will not strongly alter the conclusions of §3.2. 
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For the case of an azimuthal magnetic field, the above breaks down and a separate 
analysis is needed. There is, already at the start, a small parameter in the problem: m/k z R. 
The requirement that that resistivity be negligible is not satisfied merely by r\ B <C u, but by 
the tighter restriction 

T] B < 3u(m/k z R) 2 . (62) 
This is seen explicitly in the dispersion relation: 

m 2 \ ( 9 „ , 2 m 2 



Y + Y [VBk z + 3u—j + 7 \k 2 + Sr] B uk z ^ j + 

There can be no instability unless the constant term is negative, which requires the inequality 
(62) to hold. The asymptotic regime of validity for this condition is thus 

1 « ^ « {W) 1 ' 2 . (64) 
m 

There are further restrictions, however, if the maximum growth rate is not to depart signif- 
icantly from its value in equation (41). If we regard the dispersion relation as an equation 
for m 2 /R 2 , it may be shown that it has a well-defined solution for 7 — > 7 maa; , only if 

r) B k\ < \dtt 2 /d\nR\ . (65) 

Similarly, if we regard the dispersion relation as an equation k z , the restriction is found to 
be 

3u(m/R) 2 > \dn 2 /d\nR\ . (66) 

The last three equations may be combined to give the final form for the asymptotic domain of 
the axisymmetric dispersion relation (63) in which our solution for 7 max remains unaffected: 



k z R 
m 



{-) 


1 dn 2 


\m J 


tjb dlnR 



1/2 

< {3V) 1/2 . (67) 



4. Discussion: Origin of Galactic Magnetic Fields. 

The instability presented in this paper is generic and powerful. Its most rapidly growing 
behavior is exhibited at large vertical wavenumbers, and even in the presence of finite resis- 
tivity, it is not easily quenched. Long radial and azimuthal wavelengths are highly unstable, 
imparting large scale coherence in the disk plane already in the linear stages of instability. 
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In this section we discuss the possibility that this mechanism could be an important part of 
the process that gives rise to galactic magnetic fields. 

Generally, galactic magnetic amplification processes are divided into two principal cate- 
gories (Beck et al. 1996): (1) differential wind-up of a primordial field (Kulsrud 1986); and 
(2) dynamo amplification of a similar seed field. The latter category itself consists of at least 
two very distinct mechanisms; (2a) a classical aQ dynamo (Parker 1979), and (2b) small 
scale turbulent amplification (Schekochihin et al. 2004). The weak field limit of these mech- 
anisms are thought to be kinematic in nature. It is difficult to see how the process outlined 
in this paper could be unimportant, at least quantitatively, to any of these processes. 

Simple wind-up of a seed field by differential rotation fails to incorporate properly 
the true dynamics of either the MRI or the magneto-viscous instability described here. 
Differential rotation and essentially any field geometry is highly unstable, resulting in large 
radial motions and ultimately MHD turbulence. The shearing of radial fields of course readily 
occurs, but it is not the primary amplification process of a differentially rotating system. 

Traditional aQ galactic dynamo models rely on the presence of favorable properties of 
interstellar turbulence to generate the required mean helicity. Curiously, these theories of 
magnetic field amplification generally do not incorporate weak-field MHD instabilities, de- 
spite the latter's obvious potential benefits (e.g. Brandenburg et al. 1995; Hawley, Gammie, 
& Balbus 1996). Conversely, without something like the MRI or the magneto- viscous insta- 
bility, it is not so obvious that turbulence will generally amplify the field, at least at low to 
moderate values of V . An example of dissipative behavior is seen in a simulation of Hawley 
et al. (1996). The combination of hydrodynamical shear-layer turbulence plus a magnetic 
field lead not to dynamo activity, but to field energy loss. This result emerged despite the 
fact that calculation was done not in the kinematic limit, but fully in the MHD regime. 
However, the combination of local Coriolis, tidal, and Lorentz forces dramatically altered 
the dynamo properties of the ensuing turbulence: rapid and significant field amplification 
was observed, driven by the MRI turbulence, in a hydrodynamically stable background. 

More recently, the notion that non-helical homogeneous small scale turbulence may 
play a key role in galactic dynamos, particularly in the early stages, has been investi- 
gated in a series of numerical simulations (Schekochihin et al. 2004; see also Zeldovich 
et al. 1984, Kulsrud & Anderson 1992), which include a scalar viscosity and study the non- 
kinematical regime. The idea is that certain magnetic field configurations (termed "winning" 
by Schekochihin et al.) align themselves smoothly with the stretching direction of the strain 
tensor of the turbulent flow, but fluctuate along the corresponding null axis, so that the 
work done on the field by fluid element stretching is not undone by relaxation. A sort of 
turbulent ratchet thereby ensues, growing the field with enormous efficiency at small scales. 
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Numerical simulations of homogeneous white noise forcing conducted by Schekochihin et al. 
2004 resulted in exponential field amplification, but, interestingly, only in the regime of large 
V . What the relationship is between these small scale structures and a galactic scale field 
has yet to be established. 

Further discussion of the pros and cons of this mechanism would take us too far afield, 
but we may note that the generation of coherent magnetic field structure on scales larger 
than the disk scale height requires time scales at least as long as 1/Q. And while it may 
well be that intrinsic interstellar turbulence plays a role in the amplification of galactic 
magnetic fields, we emphasize here a phenomenon that is certainly unavoidable: differential 
rotation. In fact, there is no reason to restrict oneself to galactic radius scales: if sub-galactic, 
differentially rotating turbulent vorticies are also present, magneto- viscous modes should be 
seen on these scales as well. The key point is that a kinematic approach will miss this 
process, which is intrinsically MHD. 

One compelling origin of the seed field for the galactic amplification process is the 
stars themselves (Biermann 1950, Rees 1993), in which both battery processes and a truly 
powerful, rapid dynamo are likely to be present. A 0.1 G azimuthal surface field diluted 
in a stellar wind to interstellar scales would give rise to a seed field of ~ 2 x 10~ 9 G. With 
T ~ 10 4 K, n ~ 1 cm" 3 , one finds 

LO ci T ci ~ 10, V ~ 10 H 

and the regime of the magneto- viscous instability is valid. The ratio of gas pressure to mag- 
netic pressure is 3.5 x 10 7 , an extremely weak field by this measure. The linear amplification 
factor per orbit of the magnetic energy, as noted in §3.2, is huge: 5.2 x 10 7 . The makings of 
an MHD dynamo would seem to be present. 

If a combination of magneto-viscous and magneto-rotational processes is to be a viable 
candidate for galactic field production, it needs to be shown that (1) the bulk of the field 
energy emerges in the largest scales; and (2) the saturated field energy density can grow 
to levels comparable to the thermal pressure. Definitive answers will require a numerical 
treatment, but in the meantime the following points may be considered. 

A magnetic energy spectrum dominated by the largest scales appears to be a univer- 
sal outcome of MRI simulations (e.g. Hawley, Gammie, & Balbus 1995), despite the fact 
that the most unstable modes occur at scales much smaller than the scale height. In such 
calculations, the available dynamical range is limited. To the extent it can be measured, 
however, the inertial range is Kolmogorov-like in the magnetic energy. In hydrodynamics, 
Kolmogorov scaling for the kinetic energy power spectrum is universal, if the dissipation rate 
per unit volume is the sole constant characterizing the cascade, as often seems to be the case. 
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Universal processes may also be at work in MHD turbulence (Kraichnan 1963; Goldreich & 
Sridhar 1997). An important caveat, however, is that most MHD turbulence studies have 
been based on wave-wave interactions (which may be of secondary importance in linearly 
unstable rotating systems), and they have yet to address the role of a magnetized viscosity. 

Indeed, a magnetized viscosity may be crucial to understanding why interstellar fields are 
thermal strength (or even slightly above), whereas all numerical MRI simulations to date have 
yielded subthermal saturated field strengths. It has been argued (Balbus & Hawley 1998) 
that the outcome may be well be sensitive to V. The point is that large scale reconnection 
proceeds relatively easily in simulations when the resistive scale is comparable to or larger 
than the viscous scale. In "ideal MHD," both scales are grid based. But matters are likely 
to be very different if the viscous scale is, say, five orders of magnitude larger. In that case, 
enormous viscous stresses would occur in the course of setting up a reconnection front, and 
prevent its formation. With reconnection stifled, MHD turbulence could amplify the field to 
its natural dynamical limit: the thermal energy density. Beyond this point, buoyant effects 
would make it difficult for a suprathermal field to remain in the disk and be further amplified 
(Parker 1979). 

Large Prandtl number simulations have in fact recently begun (Schekochihin et al. 
2004), and treatments of the anisotropic Braginskii viscosity and conductivity have yet to 
be attempted. With the resistivity scale hidden below the viscous dissipation scale, there is 
no effective small scale sink for the magnetic field, and is therefore not surprising to note 
that Schekochihin et al. find that the magnetic energy spectra increases with wave number 
on subviscous scales, before it is ultimately cut-off. By way of sharp contrast, the ideal MRI 
simulations noted above find a monotonically decreasing energy spectrum for the magnetic 
field (Hawley, Gammie, & Balbus 1995). 

The magneto-viscous instability is a powerful and general mechanism to amplify weak 
magnetic fields in galaxies, or in hot dilute plasmas more generally, and is worthy of more 
detailed study. Of particular interest would be a suite of numerical simulations designed 
to isolate the differences between high Prandtl number turbulent dynamos relying on MHD 
and differential rotation, those based on random forcing only, and those containing both. 

5. Summary 

When the ion Larmor radius is less than its collisional mean free path, the form of the 
viscosity is altered. Crudely speaking, angular momentum is transported only along magnetic 
lines of force. This causes conjoined fluid elements in a differentially rotating system to 



-20- 



separate radially, dragging field lines with them. This aligns the field lines along the angular 
velocity gradient, causing an increase in the diffusive angular momentum transport, and the 
process is highly unstable. The instability has previously been studied in its plasma kinetic 
guise (Quataert et al. 2002; Sharma et al. 2003), and we find quantitative agreement in areas 
of overlap with this approach. The magneto- viscous instability does not rely directly upon 
magnetic stresses; the field serves merely to channel angular momentum transport along 
its lines of force. This is effective even at exceedingly weak magnetic field strengths. The 
maximum growth rate of the instability is given by equation (41), and exceeds the Oort-A 
value of any disk stable by the hydrodynamical Rayleigh criterion. Numerical simulation of 
the nonlinear resolution of the instability promises to be a challenging problem. Preliminary 
studies on the related magneto-thermal instability, have, in fact, begun (J. Stone 2004, 
private communication). 

The long mean free path regime is appropriate to nonradiative accretion flows around 
black holes (Quataert et al. 2002; Sharma et al. 2003), or to the interstellar medium of 
disk galaxies, our interest in this work. The magneto-viscous instability is thus a candidate 
for amplifying very small galactic seed fields into thermal strength fields. The turbulent 
MHD power spectrum, if it may be extrapolated from the moderately resolved, relatively 
low Prandtl number regime of previous examples, would evidence most of the power on 
large scales. On the other hand, the large Prandtl number regime, which characterizes 
the interstellar medium, ensures that the resistive dissipation scale is much smaller than the 
viscous scale. This may make it possible to grow thermal strength magnetic fields, despite the 
subthermal fields obtained in numerical MRI simulations, which correspond to the opposite 
Prandtl regime. But, at the same time, it is possible that the high Prandtl number regime 
may skew the power spectrum toward smaller scales (Schekochihin et al. 2004). 

A very general gasdynamical investigation of the instability, including its behavior for 
arbitrary wavenumber directions, and the dynamical stress of field line tension, will be pre- 
sented in a forthcoming paper (Islam & Balbus 2004, in preparation). 
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